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We study the problem of high temperature Bose-Einstein condensation (BEC) of atom-light po- 
laritons in a waveguide cavity appearing due to interaction of two-level atoms with (non-resonant) 
quantized optical radiation, in the strong coupling regime, in the presence of optical collisions (OCs) 
with buffer gas particles. Specifically, we propose a special biconical waveguide cavity (BWC), per- 
mitting localization and trapping of low branch (LB) polaritons imposed by the variation of the 
waveguide radius in longitudinal direction. We have shown that critical temperature of BEC oc- 
curring in the system can be high enough - few hundred Kelvins; it is connected with photon-like 
character of LB polaritons and strongly depends on waveguide cavity parameters. In the case of a 
linear trapping potential we obtain an Airy-shaped polariton condensate wave function which, when 
disturbed out of equilibrium, exhibits small amplitude oscillations with the characteristic period in 
the picosecond domain. 

PACS numbers: 67.85.Jk, 05.30.Jp, 42.82.Et, 42.50.Ct 



I. INTRODUCTION 

The investigation of quantum and statistical proper- 
ties of Bose-gases in low and especially in one dimension 
(ID) has evoked indefatigable interest in atomic optics 
and condensed matter physics for the last few decades, 
see e.g. In particular, at finite temperatures 

for a ID or 2D ideal Bose gases a true Bose-Einstein 
condensation (BEC) can only be reached in the pres- 
ence of a suitable trapping potential [||, and the criti- 
cal temperature for the phase transition depends on the 
shape of trappingpotential, which is usually harmonic 
in practice, cf. 13:] and [4]. Low-dimensional systems 
have been studied using atoms in highly deformed traps 
[Io| where effects of dimensional reduction become im- 
portant [ll|. Alternative systems are bosonic quasipar- 
ticles, where light and matter are coupled in a coher- 
ent way - see e.g. 12]. Here quantum and statistical 
properties of light Bose-quasiparticles, like excitons [l3j], 
magnons [3] and polaritons (see e.g. [L5rll8| ) have been 
considered. For example, exciton-polaritons occurring in 
quantum well structures placed in microcavities can be 
treated as a 2D gas of bosonic particles having an effec- 
tive mass which is many orders smaller than the mass of 
an electron in vacuum. This allows to study relatively 
high-temperature phase transitions in low-dimensional 
bosonic systems. Recently, the evidence of a Kosterlitz- 
Thouless phase transition - cf. (2j, superfiuid behaviour 
of exciton-polaritons in such systems, has been reported 
by many labs, [TH, [l6[. However, in the current semi- 
conductor structures the thermalization time is in the 
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picoseconds domain and comparable with the particle 
lifetime, thus we deal here with non-equilibrium "con- 
densates" [17j, for which dissipative and optical pump- 
ing effects play a crucial role, cf. 18]. The characteristic 
temperatures for these condensates are a few ten Kelvins, 
which is far above the atomic condensates, but still far 
below room temperature. 

Recently room-temperature Bose-Einstein- 
condensation of photons has been observed [H, l2fjj . 
where the photons are confined in a 2D curved-mirror 
optical microresonator filled with a dye solution. Ther- 
malization of the photon gas is established by thermal 
contact between the photons and the dye solution, using 
repeated absorption and re-emission processes in the dye 
solution [21| . Frequent collisions, on a timescale much 
faster than the excited state lifetime lead to a decoupling 
of photons and dye molecules, thus the relevant particles 
are not polaritons but photons (see [l9|, in these 
experiments. 

In the present paper we discuss a different approach to 
reach a high (room and beyond) temperature phase tran- 
sition with mixed matter-field states - polaritons in an 
atomic medium. In particular, dressed-sate polaritons, 
where a light field is coherently and strongly coupled to 
a two-level atom, leading to a bosonic quasiparticle (in a 
suitable limit) are attractive candidate due to their po- 
tentially longer lifetime, cf , [22|, [23[ . In this system ther- 
mal equilibrium of coupled atom-light (dressed) states 
can be achieved experimentally within a nanosecond do- 
main and is limited by the natural lifetime of the two- 
level atoms. Here, optical collisions (OCs) with buffer 
gas atoms can lead to thermalization, experimentally ev- 
idence for a thermal q uasi- eq uilibrium of coupled atom- 
light states has been found 0, HH] . The obtained time 
of thermalization have been about ten times shorter than 
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the natural lifetime at full optical power. Here, the lim- 
ited available power in the light field prevented a full 
thermalization. To overcome this problem, in [2(| H3] 
special metallic waveguides of various configurations with 
the length up to a few millimeters have been considered 
for trapping the polaritons inside, similar to waveguides 
and closed resonators examined for the confinement of 
microwave irradiation, cf. [281 ] . The lifetime of photon- 
like polaritons trapped in the waveguide can be longer 
than the thermalization time, and is mainly determined 
by the cavity Q-factor. Thus we expect in such a waveg- 
uide that for a large and negative atom-field detuning 
a high-temperature phase transition to a superradiant 
state of the polaritons can be reached (27| . 

In the paper we study thermodynamic and critical 
properties of dressed state polaritons trapped in a bicon- 
ical waveguide, i.e. metallic microtubes with different 
geometry 30]. In this case one can study the thermo- 
dynamics of a ID quantum gas, which exhibits a high 
temperature phase transition to a BEC in the case of 
a trapping potential which is more confining than har- 
monic Q. For this reason we consider waveguides with 
a biconical shape, which leads to a trapping potential 
in propagation direction. Throughout the paper, we as- 
sume the strong coupling limit, i.e. the eigenstates of 
the coupled atom-light system are treated as polaritons. 
In Sec. II we describe the waveguide cavity model for 
confinement of the optical mode. To be more specific, 
we consider field properties in an empty lossless metal- 
lic biconical waveguide cavity (BWC) where the photons 
are confined in transversal and trapped in longitudinal 
dimension. In Sec. Ill we give a quantum description 
of the atom-field interaction in the waveguide, and the 
problem of specifying waveguide parameters in the limit 
of weak trapping is considered. Thermodynamic prop- 
erties of the ID ideal gas of photon-like polaritons are 
examined, where we find a phase transition to a BEC 
state at experimentally feasible temperatures. In Sec. 
IV we discuss the dynamical properties of condensed LB 
polartions. In the conclusion, we summarize the results 
obtained. 



II. BICONICAL WAVEGUIDE CAVITY 

We start by describing the problem of photon trap- 
ping in the empty BWC, sketched in FigfTJ We assume 
metallic boundary conditions, and for the case of sim- 
plicity the cavity is assumed to be lossless. The radius 
of the cylindrically symmetric waveguide depends on z 
following 

with F (0) = 1 and F(|2:|) being a function increasing 
monotonically with \z\, i.e. ^ > for z > 0, similar 
to bottle resonators based on glass materials (and with 
relatively large diameter), see Ref. (26|. We here focus 




FIG. 1: (Color online) Biconical waveguide for photon (and 
polariton) trapping. Shown is a waveguide with a radius vari- 
ation following the form R(z) — Rq (1 — a \z\), which is ob- 
tained from eq. (|13|l for v = 1 in the limit of a small value of a. 
In this case the radius variation within the limited waveguide 
length is roughly linear. The radius at the waveguide center 
z — is taken to be close to A/2.61, yielding a low frequency 
cutoff close to the atomic resonance. 

on waveguides with a maximal diameter that is close to 
the wavelength of the optical field. In this limit it is 
possible to guarantee a spacing of the transverse modes 
above the thermal energy, which allows for the trans- 
verse mode quantum numbers to be frozen. The system 
becomes effectively one-dimensional, with a continuum 
of longitudinal modes above a low frequency cutoff. For 
a suitable density of states BEC then becomes possible 
in this one dimensional situation. For a related, though 
two-dimensional, situation see [13] • We consider biconi- 
cal waveguides, because fabrication of a conical hole can 
be done by e.g. ion beam etching or laser drilling, and 
then combining two of those holes to get the required 
cavity. Alternatively, self-assembling nanotubes can be 
used, which have been shown to be suitable for imple- 
mentation in the experiment [26j . 

We use the adiabatic approximation, assuming that 
R{z) is slowly varying with the z coordinate, i.e. that 
the condition 

—r 1 «! ( 2 ) 

az 

is fulfilled. In other words, we suppose that the angle 
between waveguide and the z-axis is small throughout 
the length of the waveguide - see FigfTJ 

The field properties in such a system can be described 
by the wave (scalar Helmholtz) equation for the vector- 
potential ^>(x,y^z) (we consider TM-modes of polariza- 
tion only) - cf. |32j |: 

A<Z>(x,y,z) + k 2 <Z>(x,y,z) = 0, (3) 

where we assume that the field evolves in time as e lut , 
k = lu/c is the wave vector. We can interpret $(a;, y, z) = 
<S>i(x, y, z) as a photonic wave function in the waveguide, 
where the index I numerates a set of quantum numbers 
corresponding to the solution of Eq.([3|) in terms of special 
functions. 

Both electric (E p ,E v and E z ) and magnetic 
(Hp, H v , H z ) field components can be found by using 
standard expressions for the <&(x, y, z)-function deriva- 
tives in cylindrical coordinates - see e.g. 32]. For TM- 
modes the z component of the magnetic field is H z = 0, 
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and the metallic boundary condition imposes that the 
electric field E has only a normal component on the 
waveguide surface. 

Taking into account axial symmetry of the waveguide 
in Fig Q] it is efficient to rewrite Eq. {[3j) in cylindrical co- 
ordinates z, p and if as: 



1 d ( <9$\ <9 2 $ 1 d 2 $ 
P- 



+ fc 2 $ = 0, (4) 



p dp \ dp J dz 2 p 2 dtp 2 
where the wave vector k is defined as: 

k 2 = k 2 ± + k 2 z . (5) 
We are looking for the solution of Eq. in the form: 

$(p,z,p) = *(*)e(p > z) e * n * (6) 

where m is an integer (azimuthal quantum) number, and 
the functions 9(p, z) and ^(z) comply with the equations 



d 2 <3 199 ,, 2 , , m 2 



G> 2 *(z) 

dz 2 



k 2 z M>(z) = 0. 



(7) 



(8) 



In Eq. ([7]) the terms containing derivatives on longitudinal 
coordinate z are omitted due to the adiabaticity condi- 
tion ([2]). The solution for the radial distribution 8(p, z) 
can be expressed in terms of Bessel functions, 



0(p, z) = J m (k±p), 



(9) 



where the wave vector component fcj_ is quantized in this 
case. Omitting lengthy but straightforward calculations 
for the fc_L component of the wave vector for TM-modes 
in the waveguide, we obtain: 



k±, m p{z) — ttt r — k^~ ' 



R(z) 



F(\z\) 



(10) 



where k 



(0) 



„_i_ mp — Sjj^- is a transversal wave vector compo- 
nent; g mp is p-th zero of the Bessel function J m (x) of 
the m-th order (33{. In (TTOt we have used the metal- 
lic boundary conditions for electric field components, 
Jm(k±R(z)) = 0. For p = 1 and m — we find the cavity 
transverse ground state mode, yielding k±,oi( z ) = "^jy 
with 3oi ~ 2.405. 

Equation ([5]) implies a dispersion relation of the form 



cjph = ck — c<tj ^^mp^d^Dj + k 2 . In the approxima- 
tion k±^ mp (z) 3> k z we arrive at 



W p h 



ck 



(0) 

_L,mp 



n\A) 



2k 



(0) 

_L,mp 



F(\z\) 



(11) 



In this limit the second term is much smaller than the 
first term and we assume that in the smaller, kinetic en- 
ergy term we can set -FX|z|) = 1. This approximation 



is valid in the limit of a not too large overall diameter 
variation of the waveguide. We arrive at a photon energy 



E p h = m p hc 



(hk z y 



2m 



P h 



V ph (z), 



(12) 



where we have defined an effective photon mass m p h = 
(with w cu toff = m-phC 2 /H as the low frequency cut- 



off) and an effective photon trapping potential Vph(z) 

hckjp (F(\z\) — 1) in analogy to [l9j]. The system be- 
comes formally equivalent to a one-dimensional ideal gas 
of massive particles moving along the z-axis under con- 
finement in the potential V p h.(z). 

In the following we assume a dependence 



F(\z\) = l + a\z\ l 



(13) 



corresponding to a waveguide diameter variation along 
the z-axis of R(z) = Rq/ (1 + a \z\ u ), with v ^ and 
a ^ 0, we arrive at a potential of the form V p h(z) = 

Hckjp a | z | v = m p hC 2 a \z\" . True BEC is expected to 
be possible in the ID system, if the potential is more 
confining than parabolic, i.e. v < 2 3|. We here 
mostly are interested in the case of v = 1, for which 
R{z) — Rq/ (1 + a\z\), which yields a linear trapping 
potential following V p h(z) — m p hc 2 a \z\, see Fig. [TJ For 
the case of a relatively small variation of the waveg- 
uide diameter over the length I of the waveguide, i.e. 
\al <C 1, one may use a linear variation of the diameter 
R(z) R (1 — a |z|), which gives the desired potential 
in first order approximation 26| . Such a purely biconical 
design can be easier to fabricate experimentally. 

Given the dispersion relation (|12l) the longitudinal 
wavefunction 9(z), Eq.©, follows the Schrodinger equa- 
tion 



d 2 ^{z) 2m ph 
dz 2 h 2 



(E-U ph \zn*(z) = 0, 



(14) 



where E = E p h — m^c 2 is the shifted photon energy in 
the cavity, and we made definition [7 p h ~ rn p hC 2 Q!. 

The solutions ^> n (z) of the Schrodinger equation (fl~4|) 
can be given in terms of special functions, see e.g. [33j j . 
In the quasi-classical approach the energy spectrum of 
photonic field in the waveguide can be obtained from fa- 
miliar Bohr-Sommerfeld quantization principle and reads 
(see e.g. (H) 



En 



hw v {n + 1/2)' 



(15) 



with the characteristic frequency of photonic "particle" 
oscillation, i.e. trapping frequency in the cavity 



ttU, 



l/u 

ph 



\23/ 2 rnJ 2 h( 2 -")/ 2 "I{v) ) 
Here we have defined I{v) = £ (1 - |t|") 1/2 dt. 



(16) 
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A characteristic length for the waveguide cavity can 
be obtained by looking at the turning points z = ±z c for 
the photonic ground state mode, given by 



z c = (E/u ph y /u = 



\ll(v) 



d V) (17) 



l/(u+2) 



where d v = ( ^ — 2 1 defines a characteristic lon- 

gitudinal scale of localization of the ground state (n = 0) . 
From this we can define the condition for the validity of 
the quasiclassical approach, given as (cf. [35]) 



k z du » W2) 1/3 . 



(18) 



Relation (|18|) can be satisfied for a given k z by adjusting 
the waveguide parameter a, cf. 



III. THERMODYNAMICS OF POLARITONS IN 
THE WAVEGUIDE CAVITY 

The description of atomic polaritons in a small vol- 
ume cavity can be given in the similar way as repre- 
sented in [2^,[23], i.e. by using Holstein-Primakoff trans- 
formation for atomic excitations. The Hamiltonian of 
the total system of atom and quantized field is given by 
H = -ff r ad + H at + -Hint, where H rll d characterizes non- 
interacting photons in the waveguide, H at is a Hamilto- 
nian of atomic ensemble, and Hint is responsible for the 
interaction of 7V at two-level atoms with quantized opti- 
cal field in the cavity. In momentum representation the 
Hamiltonian H may be written as 



Hk 
2N~t 



(fk+ q $l>-q<f>k(f>k> +4>{4>{,(f>k'- q fk+ q ) , 



kk' q 



(19) 



where /g (/I) is the annihilation (creation) operator for 

the photons absorbed (or emitted), </>g (</>t) is the an- 
nihilation (creation) operator that characterizes excita- 
tions (polarization) of a two-level atomic ensemble and 
obeys usual commutation relation for a Bose system, 

^ p ° b J — — ] is the collective atom-field intcr- 

action strength, p a b is atomic dipole matrix element, 
Vm is an effective volume of mode occupation within 
the region of atom-field interaction that can be defined 

as Vm = f — ffi^r ,,9\ d?r, where max ( |$(r)| 2 ] is 

Jcavity max (|*(r)| 2 ) ' V / 

maximal value of the square of the wave function, cf. 
[27l [3l| . Here we have assumed that the dipole matrix 
element <p a b is independent of the mode. This simplifica- 
tion is justified, as the exact transition rates do not influ- 
ence the thermodynamic properties of the system. The 



nonlinear part of the Hamiltonian H , i.e. the last term 
in (|19[) characterizes two-body polariton interaction pro- 
cesses due to atomic saturation effects. The polaritonic 
dispersion is determined by the photon dispersion cj p h(fc) 
and the atomic excitation dispersion, which is described 
by 



W at = W at (fc) = UJ + 



hk* 

2m„ 



(20) 



where ujq is the atomic transition frequency. 

If a quantum field intensity is not too high (the num- 
ber of photons is essentially smaller than the number of 
atoms), we can assume that the corresponding dispersion 
relation for polariton states is not modified compared to 
the uncoupled case. Thus, we can use a polariton basis 
to diagonalize the total Hamiltonian in (fT9"| by using the 
unitary transformation 



k<Pk, 



l 2,k - X k<t'k - Cfe/g, 



(21a) 
(21b) 



where the introduced annihilation operators E, g, E 2 g 
characterize polaritons in the atomic medium, corre- 
sponding to two types of elementary excitations which 
in low density limit satisfy usual boson commutation re- 
lations, cf. [23] - Parameters Xk and Ck are real Hopfield 
coefficients satisfying condition X\ + C\ — 1, which de- 
termines the contribution of the photon (Ck) and atomic 
excitation (Xk) fraction to the polariton annihilation op- 
erators (|2Tj) according to 



Xh — 



1 



1/2 



y/2 \ V 4 ^ 2 + S k 



C k = -= 1 



Sk 



V2 I V 4k2 + 5 l 



where Sk = w ph - w at « A 
quency mismatch. Here A = 



fikz 



1/2 



V ph {z) 



(22a) 
(22b) 

2^ t- —5— is the fre- 
wl — wo is the atom-field 
detuning, where lol is laser light frequency which is taken 
to be close to w cuto ff • In the presence of photon trapping 
Xk and Ck parameters depend on waveguide longitudinal 
coordinate z. 

We confine our analysis to polaritons of the lower 
branch (LB polaritons), and we ignore the effects of in- 
teractions between lower and upper polariton branches. 
Taking into account quasiclassical approach (cf. (|18)l) for 
the photonic field in the BWC we find the general con- 
ditions for the observation of a BEC of (lower branch) 
polaritons, 



AE n < k B T<^hfl R0 , 

iV2 



(23) 



where fi^o = (A 2 + 4k 2 ) ' " is zero-momentum Rabi 
splitting frequency between the lower and upper polari- 
ton branch. The first constraint in (|23[) represents the 
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condition for a quasiclassical limit where the energy spac- 
ing AE n ~ fiiUy (see (|15|l) of the quantized photonic 
states is essentially smaller than the thermal energy and 
the states may be treated as a continuum, cf. where 
the energy levels are populated according to a Bose- 
Einstcin distribution in thermal equilibrium. The second 
condition implies that the thermal energy is not enough 
to excite lower branch polaritons to the upper branch, 
essentially allowing us to neglect the upper branch. 

Taking Eqs. ([2"T]). ([25]! into account from Eq.JTOJl we 
can obtain LB Hamiltonian i?LB in the form 

#LB = h ^ n k E l k E 2X + Yl U kk'q El 2,k+ q E 2,k'-q E ^^2M', 

k kk 'q 

(24) 



where f2r = j [u! at + <^ph — Qr] determines the disper- 

1/2 

sion relation for LB polaritons; D,r = (8 2 + 4k 2 ) is 
the Rabi splitting frequency that determines the gap 
between upper and lower states. The gap is minimal 
and equal to Qro taken for k z — at the center of the 
trap (waveguide) at z = 0. In Eq. (|2~4"|) we give the defi- 

niti0n U kk'g = { C \k + q\ X k> + C 'k' X \k + q\) X \k'-q\ X k 

that determines two-body polariton-polariton scattering 
processes. The effective mass of the LB polaritons is 
found to be 



dk 2 



2m at fflphO & (z) 



fe,=o/ (m at -I- TOp h ) Vt Rz (z) - (m at 

I 



m ph ) (A + V ph (z) H 



(25) 



where Qr z (z) — 



4k' 



1/2 



the 



dependent Rabi splitting frequency. 

To be more specific we examine the interaction be- 
tween a quantized field and rubidium atoms, which are 
treated as a two-level system. The transition frequency 
is taken as the weighted mean of the rubidium D-lines, 
loq/2it ~ 382 THz, and the laser is taken to be red de- 
tuned by |A|/2tt ^ 11 THz , A < In the 
experiments reported the thermal energy (ksT) for the 
atoms at ambient temperatures (T = 530 K) is about 
the frequency of Rabi splitting f!,Ro/2Tr ~ 11 THz in en- 
ergy units, thus condition (|23|) has not been completely 
achieved yet [24Tj27j . We consider the perturbative limit 
when fl , k < |A|, that is fl RQ /2TT sa |A| /2tt ^ 11 THz. 
For negative detuning (A < 0) the LB polaritons with 
relatively small momentum k z at the center of the trap 
are photon- like with the mass m po i « m phf7~fra' — 
2.8 x 10" 36 kg, that is S k ^ A, X k w 0, C k 1, and 

The role of polariton-polariton scattering processes in 
Eq. (f2"4")) is negligibly small in this limit. At the bottom of 
the dispersion curve the polariton scattering parameter 



^kk'q 



behaves as Un 



Hence in the 



limit of a large atom-light detuning |A| (ratio k/ |A| is 
about 0.057, cf. [25|, [27|) and for a macroscopically large 
number of atoms (N a t ^> 1) the gas of LB polaritons can 
be treated as an ideal one. 

The photonic fraction of the LB polaritons - Hopfield 



coefficient C, 



2 — 

— 



r 2 



and the mass of LB polari- 



tons TOp i taken in the limit of zero momentum k z — 
as a function of z are presented in Figf2^,,b. With in- 



creasing coordinate z, LB polaritons become more atom- 
like, and present half matter - half photon quasi-particlcs 
(X (z) = Cq(z) = 2~ 1 / 2 ) with the mass m po i « 2m p h at 

the distance Zi/ 2 — (^|A| /U- P \ i ) 1 ^ v , where photonic and 
atomic dispersion lines cross (see Fig|4ji). From definition 
(f2l))) it follows that the mass of polaritons becomes z in- 
dependent if condition V p h (z) <C ?i|A| is satisfied. Since 
Vph {d v ) ~ fiw„, the inequalities (|23p . already involve the 
above condition. In this limit an effective width of the 
photonic mode localization d v is essentially smaller than 
the characteristic length £1/2- The ratio Z]y 2 /d„ for dif- 
ferent trapping power is outlined in the inset in Fig|2]3. 
For v — 1 the region of photon mode localization is more 
than ten times shorter than characteristic length 2^/2, 
allowing to treat the LB polaritons as nearly photonic. 
The variation of the polariton mass, see Figj2]D, is small 
for z <C Zxiii i.e. under condition (f2"3"|). allowing to treat 
the LB polaritons as particles with constant mass. 



Figure [3] shows the characteristic trapping frequency 
uj v as a function of parameter v for the cavity, where red 
(dashed) curve corresponds to waveguide a-parameter 
used in Fig. [5] For v > 1 the condition (|2"3")) can be ful- 
filled for any reasonable value of a-parameter by choos- 
ing the atom- field detuning |A|. On the other hand, uj v 
increases within the domain of < v < 1 and the fulfill- 
ment of condition (|23[) mainly depends on the given value 
of curvature of the waveguide radius, i.e. on a-parameter 
and on experimentally accessible atom-field detuning | A|. 

The results for the Rubidium system allow us to fur- 
ther simplify the Hamiltonian for the LB polaritons: At 
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FIG. 2: (Color online), (a) - The polariton photon frac- 
tion (coefficient Cq) and, (b) - normalized LB polariton 
mass m P oi as a function of z for trapping power parameter 
v = 0.5, 1, 1.5. The inset shows dependence of z 1 i 2 /d u ver- 
sus parameter v. The parameters are A/2n = —11 THz, 
Ro ~ A/2.61 « 0.3 /zm, a 17 " = 0.0005 /urn -1 . 




FIG. 3: (Color online). Dependences of photonic trapping 
frequency uj v /2n versus trapping parameter v. Blue solid line 
corresponds to a 1/V = 0.005 pm' 1 , red dashed - to = 
0.0005 firn' 1 , green dotted - to a 1/lJ = 0.00005 pm' 1 . In the 
inset the dependence of cj„ versus parameter a 1 ^ for v = 0.5 
(blue solid line), v = 1 (red dashed line) and f = 1.5 (green 
dotted line) is depicted. 



small momenta 



<C ^_ro it approaches 



ill 



LB 



E 



21.2 



2m pol 



(26) 



where J7(z) is an effective trapping potential for polari- 
tons defined as 



U(z) 



(HA + V ph (z)) 2 + Ah 2 n 2 + Ml 



no 



(27) 

The last term in brackets of Eq. (|27|) specifies the mini- 
mal level of potential energy U (z), which is equal to zero, 
U(z)\ z=Q = at the center of the trap. From Eq. (|2"Tf un- 
der the condition (l23l) for photon-like polaritons one can 
obtain a simple expression for power-law potential 



U{z) 



pol 



(28) 



where U po i = m ph c 2 a (tt R o + \A\)/2Q m . 

Let us examine statistical properties of LB polaritons 
trapped in the waveguide. From the above discussion, we 
can treat the polaritons as one-dimensional ideal bosons 
confined in the potential U{z). In the quasiclassical ap- 
proximation (|23p the density of states approaches 



V2 



1(e) 



e-U(z) 



1/2 



dz, 



(29) 



l/u 



where I (e) = (jf~[) ^ s a characteristic length of lo- 
calization for LB polaritons with energy e - see Fig|4ji. 
The density of states p (e) as a function of the normal- 
ized polariton energy e/fi,|A| is plotted in FigHb for the 
approximate (dashed) and exact (solid) potential. Differ- 
ences occur at polariton energies of e ~ ft|A|, i.e. for LB 
polaritons weakly confined inside the region z < Z\/2 one 
can use power-law approximation (|28[) of the trapping 
potential U(z). 

The total number of LB polaritons N po \ is given by 



N 



pol 



N Q 



p (e) de 



exp[(e- p)/k B T}-l' 



(30) 



where No is the number of ground state polaritons, p is 
chemical potential. 

We find the critical temperature Tc, for which the 
ground state occupation becomes macroscopically by 
solving Eq. ([5H|) at p — 0. For the onset of Bose-Einstein 
condensation one finds (cf. 3) 



knT( 



B±C 



where F(v) 



^/2rri~iF(iy)T(x)C(x) 
i t 1 '"- Ut 



2v/(2+v) 



(31) 



J Q J t dt ; T(x) and ((x) are the gamma 
and the Riemann (zeta) functions taken at x = l/v+1/2, 
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FIG. 4: (Color online) (a) shows the energy of polariton trap- 
ping versus coordinate z taken for v = 1; (b) shows the de- 
pendence of density of states p (e) on normalized polariton 
energy e/h\A\. The solid curves correspond to exact trapping 
potential (|27(l . and the dashed ones are relevant to power- law 
potential (I28|l . The black solid lines correspond to atomic 
Wat and photonic oj p h frequencies. The shaded (yellow) area 
corresponds to LB polariton ground state wave function - see 
Sec. IV. 




200 400 600 800 1000 



N pol 

FIG. 5: (Color online). Critical temperature Tc versus num- 
ber of polaritons N po \. 



temperature Tc as a function of LB polariton number 
./Vp i is shown. High temperature BEC, as it follows from 
Eq. (l31[) and Figj5]can be achieved for an experimentally 
feasible number of polaritons. Since the function is 
a diverging at x = 2 (see Eq. (l3lj) and [11]), the critical 
temperature Tc vanishes for increasing trapping power 
parameter v - see Fig(5] From Eq. ([3"Tj) we find a critical 
temperature which drastically increases with increasing 
Q-parameter. However, the values of the parameter are 
limited by inequality @ for our problem. 

Notably, the thermal de Broglie wavelength At = 

r 9 fi 2 i 1/2 

m ifc B T a * ^ ne experimentally explored tempera- 
tures of the atomic gas T = 530 K is macroscopically 
large, i.e. A = 1.89 fim and comparable with the mag- 
nitude of characteristic length d\ of photonic field local- 
ization. 



IV. POLARITON BEC PROPERTIES AT v = 1 



respectively. Below the critical temperature the occupa- 
tion of the ground state is then determined by 



pol 



1 



T 



(32) 



In the experiment, the average number of LB polari- 
tons Ap ol = J^k ~ ^P h can ^ e estimated us- 
ing the photon-like character of polaritons in the pertur- 
bation limit. Notice that a polaritonic model is valid 
under the low excitation density limit for which the 
average number of photons Af p h is essentially smaller 
than the average number of atoms JV a t. Using exper- 
imentally accessible rubidium atom densities, that are 
"at = N at /V M = 10 16 cm~ 3 , cf. [24], and taking the oc- 
cupation volume Vm of the lowest photonic mode (TMoi 
mode) to be Vm ~ 0.5 /im 3 we find an average number of 
A^at = 5000 atoms in the BWC, thus limiting the exper- 
iment to a few hundred polaritons. In FigjS] the critical 



Let us examine LB polariton condensate properties at 
sufficiently "low" temperatures such as T <C Tc which is 
relevant to the linear trapping potential U(z) obtained at 
v = 1. In this case the parameter U po \ — U (z) / \ z\ may 
be physically interpreted as a force acting on polaritons 
in the waveguide cavity. Macroscopic LB polariton BEC 
properties can be found with the help of a quantum field 
theory approach. In particular, the Lagrangian density L 
for the system described by Hamiltonian (|26p looks like 



L ^ 
~ ~2 



^-df-^ at 



2 in 



pol 



dz 



where tp(z,t) = 
densate ground 



° \ip(z, t)\ 2 dz 



U pol \z\W, 
(33) 

S(z, t)^ is classical polariton con- 
state wave function normalized as 
= 1 . Lagrangian density L implies non- 



stationary Schrodinger equation for polaritons in coordi- 
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note representation in the form 



ih 



dt 



2m po \ dz 2 



+ U po i\z\ 



1>{z,t). (34) 



To obtain a stationary solution we make a substitution 
ip(z,t) = e^ lfJlt ^ h ^p(z), where /x is a chemical potential. 
Then we have 



d 2 i\){z) 



2m po i dz 2 



+ (M - f/poi \z\) i>{z) = 0. (35) 



The LB polariton ground state wave function is ex- 
pressed through the Airy function as, cf. 3(i : 



where d\ 



pol 



1.308 



Ai 



1 = 1 



(36) 



1/3 



specifies a characteristic 



2m pol £/ pol 

scale of spatial (longitudinal) localization of the polariton 
condensate for a linear trapping potential; a[ ss —1.0188 
is the first zero of the Airy function derivative. It is 
interesting to note that the characteristic length di. po i = 
di is the same for the photons and photon-like polaritons 
trapped in the BWC due to relation m po iE/poi = TO p h£/ph> 
that is true under the condition (|2"3")1 . see also Eqs. (TTTl) . 
(l28f . Since LB polaritons are completely photon- like it is 
possible to conclude that the photon mode n = reflects 
the behaviour of polaritonic probability density \ip n \ = 
l^l . The chemical potential [i can be easily found from 
Eqs. (HI])-® and given by 



fi = 1.0188 U pol di. 



(37) 



To investigate the reaction of the BEC to a small dis- 
turbance induced from the outside we study the dynamics 
of the BEC when the initial wave function has been com- 
pressed or stretched compared to the equilibrium wave 
function. In this case one expects slow oscillations of the 
BEC around the equilibrium state. The dynamics of the 
polariton BEC induced by a disturbance can be found by 
means of the variational approach for solving Eq. (|34p . In 
particular, we take the Airy trial function for the ground 
state (cf. [37[): 



tl>(z,t) 



a ib(t)\z\ 



(38) 



N is a normalization constant. In Eq. (|38|) the time- 
dependent function D{t) specifies the width of the wave 
function, and b(t) characterizes a related wave function 
curvature. Inserting Eq. ([38j) into Eq.([33|) it is possible to 

+00 

obtain an effective Lagrangian L = J Ldz by averaging 



Lagrangian density L as 



2ha[ db 



D 



h 2 a[ 1 



3 dt 6m po i D 2 2m po i 



h 2 ^ 2 _ 2U po ia' 1 



(39) 




FIG. 6: (Color online) Normalized wave function width w 
versus time t. The parameters are the following: v = 1, 
Ro w 0.3 /xm and a 1/v = 0.0005 fim . An initial condition 
is w (t — 0) = 0.75 for solid (blue) curve, w (r = 0) = 0.8 for 
dashed (red) curve, and w (r = 0) = 0.9 for dotted (green) 
curve. In all cases w (r = 0) =0. The horizontal line w = 1 
corresponds to the stationary solution of Eq. (l42[) . 



The effective Lagrangian L leads to a Newton-like 
equation for condensate width function 



d 2 D 
dt 2 



3ft 2 



3£7, 



p ol 



4 \a\ | m 2 pol D3 



ai rn 



pol 



(40) 



For further processing it is useful to introduce a new di- 
mensionless variable for the wave function width w = 
D/di and the rescaled time r = u) po it, where w po i is a 
characteristic frequency of polariton trapping defined as 
(cf. {16])) 



Wpol 



V3ft 



3 3 / 2 rj p 2 0l 



1/3 



I a'i |Ti po id 2 



2 la' 



,3/2 



ftTO 



(41) 



pol / 



The dimensionless equation for the new variables takes 
the form 



d 2 w 
W 2 



1 



- 1. 



(42) 



The solution of Eq. (|42)) in the absence of polariton trap- 
ping is w = vl + r 2 ; describing the spreading of a free 
polariton wave packet with Airy shape - cf. [36j . In 
the presence of polariton trapping the analytical solution 
of Eq. (|42p is much more complicated. Figure [B] demon- 
strates the temporal dynamics of the normalized polari- 
ton condensate wave function width w, using the Ru- 
bidium parameters described in section [TTTl (cf. [2~I1 - [26T | V 
which shows oscillations around the equilibrium solution. 
The characteristic frequency estimated for experimental 
conditions is cj po i/27r = 0.817 THz. It is important that 
the value of w po i should satisfy the quasiclassical con- 
dition (f2"3"]) discussed above. Completely neglecting the 
decay rate, the polaritonic system exhibits periodical be- 
havior similar to an atomic BEC in a harmonic trap - cf. 
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[37| . The frequency f2 osc for small amplitude oscillations 
(see Fig (5]) can be obtained by linearizing Eq. (|4"2")) around 
the stationary state w = 1 is found to be f2 osc = V^poi- 
The polariton wave function is stable everywhere in this 
case. 

V. CONCLUSIONS 

Let us briefly summarize the results obtained. We have 
considered the problem the thermodynamics for polari- 
tons emerging due to the interaction of two-level rubid- 
ium atoms with an optical field in the presence of OCs 
with buffer gas particles. We assume the optical field to 
be trapped inside a biconical waveguide filled with Ru- 
bidium atoms, and we find conditions for effective trap- 
ping of the coupled atom-light system. In particular, we 
find a phase transition to a BEC of polaritons if relation 
([23]) is fulfilled, a condition that substantially depends 
on a cavity mode structure and a-parameter characteriz- 
ing slow changing of waveguide radius along longitudinal 
coordinate. We analyze the problem of BEC formation 
for lower branch atomic polaritons trapped in the BWC 
under the quasiclassical approach (fT8|) . ([23]). and we find 
that atomic polaritons formed in the waveguide may be 
treated as ID ideal gas of bosonic quasiparticles in the 
framework of current experimental results achieved for 
the thermalization of coupled atom-light states in the 
presence of strong atom-field coupling regime. 

Even though the character of the dressed-state polari- 
tons discussed in the paper is nearly photon-like, they 
represent mixed states of photons and excitations of two- 
level atomic system. In this sense they are physically 
close to exciton-polaritons obtained in semiconductor mi- 
crocavities, cf. (15l [l6j. However, our estimations show 
that the critical temperature of the BEC phase transition 
can exceed the temperature of 530 K currently used in 



experiments, thus being high enough to be observed us- 
ing appropriate waveguide parameters and cavity modes. 
Such high critical temperatures can not be achieved for 
exciton-polaritons in the current narrow-band semicon- 
ductor microstructures due to exciton ionization effects. 
The main reason for the high transition temperature is 
the photon-like character of the polaritons, i.e. their low 
effective mass. This, together with the thermalization 
due to optical collisions, may allow to experimentally 
observe a high temperature phase transition of dressed 
state polaritons using realistic parameters. 

In a more detailed study we have investigated the case 
of a biconical cavity leading to a linear trapping potential 
of polariton gas. We are aware that the transverse deco- 
herence rate needs to be studied in more detail for the 
considered system to ensure the strong coupling limit. 
Using variational techniques we find the ground state 
wavefunction to be Airy-shaped, where the width w ex- 
hibits small amplitude oscillations in time around a sta- 
tionary state when the BEC is excited, with a character- 
istic period in the picosecond domain for current exper- 
iments. The results obtained indicate that the observa- 
tion of a corresponding thermodynamic phase transition 
is possible already in current experimental setups, given 
the described waveguide structures are prepared with the 
required accuracy. 
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